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Power Spectral Analysis of Heart Rate and 
Arterial Pressure Variabilities as a Marker of 
Sympatho- Vagal Interaction in Man and 
Conscious Dog 

Massimo Pagani, Federico Lombardi, Stefano Guzzetti, Ornella Rimoldi, 
Raffaello Furlan, Paolo Pizzinelli, Giulia Sandrone, Gabriella Malfatto, 
Simonetta Dell'Orto, Emanuela Piccaluga, Maurizio Turiel, Giuseppe Baselli, 
Sergio Cerutti, and Alberto Malliani 



In 57 normal subjects (age 20-60 years), we analyzed the spontaneous beat-to-beat oscillation in R-R 
interval during control recumbent position, 90° upright tilt, controlled respiration (n = 16) and acute 
(n - 10) and chronic (n = 12) ^-adrenergic receptor blockade. Automatic computer analysis pro- 
vided the autoregressive power spectral density, as well as the number and relative power of the 
individual components. The power spectral density of R-R interval variability contained two major 
components in power, a high frequency at -0.25 Hz and a low frequency at -0.1 Hz, with a 
normalized low frequency : high frequency ratio of 3.6 ± 0.7. With tilt, the low-frequency component 
became largely predominant (90 ± 1 %) with a low frequency ; high frequency ratio of 21 ± 4. Acute £- 
adrenergic receptor blockade (0.2 mg/kg IV propranolol) increased variance at rest and markedly 
blunted the increase in low frequency and low frequency : high frequency ratio induced by tilt 
Chronic ^-adrenergic receptor blockade (0.6 mg/kg p.o. propranolol, t.i.d.), in addition, reduced low 
frequency and increased high frequency at rest, while limiting the low frequency : high frequency ratio 
increase produced by tilt. Controlled respiration produced at rest a marked increase in the high- 
frequency component, with a reduction of the low-frequency component and of the low frequency: 
high frequency ratio (0.7 ± 0.1); during tilt, the increase in the low frequency : high frequency ratio 
(8.3 ± 1.6) was significantly smaller. In seven additional subjects in whom direct high-fidelity arterial 
pressure was recorded, simultaneous R-R interval and arterial pressure variabilities were examined 
at rest and during tilt. Also, the power spectral density of arterial pressure variability contained two 
major components, with a relative low frequency : high frequency ratio at rest of 2.8 ± 0.7, which 
became 17 ± 5 with tilt. These power spectral density components were numerically similar to those 
observed in R-R variability. Thus, invasive and noninvasive studies provided similar results. More 
direct information on the role of cardiac sympathetic nerves on R-R and arterial pressure variabilities 
was derived from a group of experiments in conscious dogs before and after bilateral stellectomy. 
Under control conditions, high frequency was predominant and low frequency was very small or 
absent, owing to a predominant vagal tone. During a 9% decrease in arterial pressure obtained with 
IV nitroglycerin, there was a marked increase in low frequency, as a result of reflex sympathetic 
activation. Bilateral stellectomy prevented this low-frequency increase in R-R but not in arterial 
pressure autospectra, indicating that sympathetic nerves to the heart are instrumental in the genesis of 
low-frequency oscillations in R-R interval. {Circulation Research 1986;59:178-193) 



THE whole history of brain electrophysiology 
has proved that rhythms can often be markers 
of normal functional states such as wakeful- 
ness or sleep, or of abnormal states such as epilepsy. In 
the present context, the general hypothesis is that 
rhythmical beat-to-beat oscillations of one cardiovas- 
cular controlled variable might provide some criteria to 
interpret the complex interplay among the neural regu- 
latory outflows. On the other hand, the existence, un- 
der stable conditions, of rhythmic fluctuations in car- 
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diovascular variables such as heart rate or arterial 
pressure has been recognized for a long time. 1 - 7 In 
recent years, the possibility of quantifying these oscil- 
lations by using computer techniques, particularly on 
the variability of electrocardiographic R-R interval, 
has aroused a growing interest. 8 " 15 

Indeed, as instantaneous heart rate depends on the 
interaction between sympathetic and parasympathetic 
efferent activities and pacemaker properties, it has 
been suggested that this analysis could lead to a nonin- 
vasive assessment of the "tonic" autonomic regulation 
of heart frequency. 12 The nonrandom components of 
R-R variability usually are assessed with spectral tech- 
niques based on the fast Fourier transform, in spite of 
the consideration that the heart rate variability signal is 
not strictly periodical, as requested by the determinis- 
tic nature of the algorithm. 16 
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In conscious dogs and in man, a high-frequency 
component (—0.25 Hz) has been consistently found in 
the power spectrum and has been interpreted as a quan- 
titative assessment of respiratory arrhythmia.' 5 One or 
two low-frequency components have also been de- 
scribed, respectively, about 0.1 Hz and 0.03 Hz. The 
former of these, the so-called 10-second period 
rhythm, has been considered particularly interesting as 
it has the same frequency of the better known Mayer 5 
waves.' 7 

As to the neural mechanisms underlying these fluc- 
tuations, vagal efferent activity has been interpreted as 
being primarily responsible for the high-frequency 
component of heart rate variability on the basis of 
experiments on vagotomized decerebrated cats, 10 con- 
scious dogs, 12 and man with muscarinic receptor 
blockade. 15 

Both vagal and sympathetic outflows were consid- 
ered to determine the low-frequency components. 1215 

The aim of this study on normal human subjects and 
conscious dogs was to assess the relative role of vagal 
and sympathetic activities in determining the variabil- 
ity in heart rate and arterial pressure at rest and during 
induced changes of autonomic regulation. A totally 
automatic computation of the autospectra and of their 
individual components was implemented, using an 
autoregressive modeling algorithm, in order to take 
advantage of its inherent better definition, compared to 
the more current spectral techniques. 18,19 

Subjects and Methods 

Fifty-seven subjects without any clinically evident 
disease (20-60 years old) were used for this study. 
They were healthy volunteers, randomly selected from 
hospital staff, medical students, and their relatives. 
Seven additional subjects consulted our hypertension 
clinic and were eventually found to be normotensive in 
the hospital environment. These subjects underwent 
24-hour continuous monitoring of electrocardiogram 
(ECG) and high-fidelity intraarterial systemic pressure 
recording 20 as part of their clinical workup for suspect- 
ed hypertension. 

All subjects included in this study had no signs of 
organ or systemic disease. Subjects smoking more 
than five cigarettes a day were excluded from the 
study. No one was taking any medication, and each 
subject was instructed to avoid beverages containing 
alcohol or caffeine after 10 pm of the day preceding the 
study and to come to the laboratory after a light break- 
fast. Studies were performed between 10:30 am and 
12:30 pm. All subjects were carefully instructed about 
the study, and all gave their informed consent. 

Tilting 

Subjects were placed on an electrically driven tilt 
table (provided with a mattress and a footrest) and 
connected to a conventional electrocardiographic am- 
plifier (lead II, Cardiol ine) and FM tape recorder (Ra- 
cal). After 15-30 minutes, allowed for stabilization, 
the ECG was continuously recorded for 30 minutes at 
rest, followed by 20 minutes after the subjects had 



been passively moved to an upright 90° position by 
electrically rotating the table. During tilting, the sub- 
jects rested on their feet, although they had been loose- 
ly strapped to the table to avoid accidental falls. 

The arterial blood pressure was recorded at the be- 
ginning and at the end of "rest" and "tilt" periods using 
a conventional sphygmomanometer. 

Reproducibility of the results was assessed by re- 
peating the study in 10 subjects after a period of 1 week 
to 12 months. Additionally, in 4 subjects, a third study 
was performed after an additional 3-4 months. 

The effects of ^-adrenergic receptor blockade was 
assessed in two groups of experiments. In 10 subjects, 
after the baseline studies, the study was repeated on a 
later day after IV bolus of 0.2 mg/kg propranolol (In- 
deral, ICI) in order to block acutely /^-adrenergic re- 
ceptor-mediated effects (Table 1). 

In 1 2 subjects, the study was repeated after 6 days of 
0.6 mg/kg p.o. propranolol (Inderal, ICI) t.i.d. in or- 
der to block chronically ^-adrenergic receptor mediat- 
ed effects (Table 1). 

Controlled Respiration 

The effects of controlled respiration were tested in 
16 subjects both at rest and during tilt. These patients 
were studied first while breathing spontaneously, and 
on a separate day while breathing at 20 acts/min, fol- 
lowing a metronome. Respiratory movements were 
monitored with an impedance respirometer (Cardio- 
line) and recorded on FM tape. Seven of these patients 
breathed through a mouthpiece connected to a spiro- 
meter (Jaeger) in order to quantify their tidal volume. 
All of these patients had been well acquainted with the 
laboratory through 1-3 training sessions, during which 
they learned to breathe quietly and to adjust their tidal 
volume to the increased respiratory rate. During spon- 
taneous and controlled breathing, transcutaneous-Pco 2 
and Po 2 (Sensor Medics) were monitored in 5 patients. 

Systemic Arterial Pressure (Invasive Studies) 

In 7 additional subjects, a microminiature Millar tip 
pressure transducer (03 French) was inserted percuta- 
neously into the radial artery of the nondominant arm 
with a Seldinger technique. 

Calibrations of the transducer, both in absolute val- 
ues of millimeters of mercury and millivolts, were 
performed just before the insertion of the catheter and 
at the end of the recording procedure. After suitable 
amplification, arterial blood pressure and ECG were 
continuously recorded at rest and during tilt on FM 
tape. These studies were the last part of a 24-hour 
protocol for suspected hypertension, based on a con- 
tinuous ambulatory recording of direct high-fidelity 
arterial pressure and ECG. 20 

Animal Experiments 

Twelve mongrel dogs (25-35 kg body weight) were 
used for this part of the study. They were lightly anes- 
thetized with thiopental sodium (10 mg/kg, IV) (Far- 
motal, Farmitalia) and, after the skin had been infiltrat- 
ed with 2% xylocaine (Byk Gulden), a cut-down was 
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Table 1. 



Effects of Age and fi- Adrenergic Blockade on the Response of R-R Interval Variability to Tilt 



Arterial pressure 

Diastol- 
Systolic ic (mm 
(mm Hg) Hg) 



R-R 
interval 
(msec) 



R-R 
variance 
(msec 2 ) 



Low-frequency 
component 

Normalized Frequency 
power (Hz Eq) 



High-frequency 
component 

Normalized Frequency 

power (Hz Eq) LF: HF ratio 



A. Rest 

20-30 yr (n = 30) 

30-45 yr (n = 10) 

45-60 yr (n = 17) 

Acute /^-adrenergic 
blockade, 20-30 yr 
(n = 10) 

Chronic ^-adrenergic 
blockade, 20-30 yr 
(i = 12) 

B. Tilt 
20-30 yr 

30-45 yr 

45-60 yr 

Acute j3-adrenergic 
-blockade, 20-30 yr 

Chronic 0- adrenergic 
blockade, 20-30 yr 



116 

±2 

116 
±4 

118 
±3 

103 
±3t 



73 
±2 

79 
±3 

78 
±2 

70 
±3 



834 
±34 

931 

±39 

926 
±27 

1076 
±26t 



4097 
±361' 

2581 
±356' 

1354 
±205 

8255 
± 1676t 



58.2 
±3.3 

62.2 
±6.3 

49.9 
±4.9 

50.9 
±5.0 



0.11 
tO.Ol 

0.10 
tO.Ol 

0.09 
tO.Ol 

0.12 
tO.Ol 



24.0 
±2.0 

26.3 
±4.3 

32.0 
±4.6 

32.7 
±4.3 



0.24 
t0.02 

0.27 
t0.02 

0.27 
tO.Ol 

0.26 
t0.02 



3.62 
±0.70 

3.69 
± 1.23 

2.58 
±0.61 

1.74 
±0.31 



106 
±2t 


66 
±2t 


1115 

±59t 


7180 
±1712t 


39.4 
±5.7t 


0.12 
±0.01 


53.4 
±5.4| 


0.26 
±0.02 


1.07 
±0.35t 


111 
±1 


79 

±2 


.687 
± 14 


t 


2642 
±217 


89.7 
t ±1.4 


$ 


0.09 
±0.01 


7.5 
±0.9 




0.29 
±0.03 


20.79 
±3.68* 


109 
±4 


78 
±3 


719 
±44' 


t 


2160 | 
±253«j' 


83.7 
±4.6' 


t 


0.09 
±0.01 


U.4 
±2.8' 


t 


0.29 
±0.06 


. 17.30 
±8.19* 


113 

±3 


77 
±2 


767 
±24 


t 


1212 
±2131 


75.7 
±4.6 


$ 


0.08 
±0.01 


16.2 
±3.4 


i 


0.23 
±0.02 


12.61 
±3.21* 


100 
±4t 


72 
±4 


929 
±24t* 


3009 
±361* 


74.1 
±7.0tt 


0.09 
±0.01 


13.0 
±2.5t 


0.22 
±0.02 


7.48 
± 1.77ft 


97 

±2ti 


69 
±4t 


951 

±35tt 


4067 
±664t$ 


66.7 
±5.2tt 


0. 11 

±0.01 


24.0 
±4.4t* 


0.22 
±0.02 


3.78 
±0.93tt 



t Value during ^adrenergic blockade significantly different (p < 0.05) from value obtained in unblocked conditions in the young age group 
TValue during tilt significantly different from value at rest (p < 0.05). ^^uup. 



performed in the region of the femoral artery. A cath- 
eter was introduced into it, advanced to the aorta, and 
secured with a suture. The wound was closed, and the 
catheter was exteriorized at the base of the neck. 

After baseline studies, 5 dogs underwent bilateral 
stellectomy. By sterile techniques, they were anesthe- 
tized with thiopental sodium (30 mg/kg, IV) followed 
by a continuous infusion of fentanyl citrate (6.2 jag/kg, 
IV) and droperidol (0.2 mg/kg, IV) (Leptofen, Carlo 
Erba). The animals were paralyzed with small (0.1 
mg/kg, IV) intermittent doses of succinylcholine 
(Wellcome) and artificially ventilated with a positive- 
pressure pump (Harvard). A thoracotomy was per- 
formed in the fourth left intercostal space, the left 
stellate ganglion and its branches were excised, and the 
thoracotomy was repaired. Seven to 15 days later, 
under similar anesthesia and through a right thoracot- 
omy in the fourth intercostal space, the right stellate 
ganglion was excised as well. 

Aortic pressure was measured with the implanted 
catheter using a pressure transducer (Statham Instru- 
ments). Aortic mean pressure was obtained with an RC 
filter with a 2-second time constant. The electrocardio- 
gram (lead II) was obtained with subcutaneous silver 
electrodes and an AC amplifier. Heart rate was moni- 
tored continuously with a cardiotachometer triggered 
by the R-wave. 



Respiratory movements were monitored with a 
pneumatic belt connected to a pressure transducer 
(Statham Instruments). Data were recorded on a multi- 
channel FM tape recorder (Racal Store 7) and played 
back on a direct-writing recorder (Brush Gould). 

Experiments were performed after a postoperative 
period long enough to allow complete recovery of the 
dogs from the operation, as judged by their normal 
behavior, body temperature, and hematocrit. Five to 7 
days usually were sufficient after the implantation of 
the pressure catheter, while 10-15 days were given 
after the thoracotomy. While the trained dogs were 
lying quietly on the recording table, aortic blood pres- 
sure, ECG (lead II), and respiratory movements were 
recorded continuously for 20-30 minutes of control, 
followed by an infusion of nitroglycerin (32 /-eg/kg per 
min, IV) for 15-20 minutes, in order to excite sympa- 
thetic activity. 21 

Data Analysis 

Off-line analysis was performed on DEC MNC 
11/23 and PDP 11/24 minicomputers. ECG, respira- 
tion, and pressure data were played back from the FM 
tape and digitized at 300 samples/sec per channel. The 
principles of the software for data acquisition and anal- 
ysis have been described previously. I3,2 ° i22 

Briefly, stationary sections of data both at rest and 
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Figure 1. Schematic outline of the computer analysis of R-R variability. From the surface ECG (top trace), the series of R-R 
intervals is calculated as a function of the beat number. This gives rise to the tachogram and, from it, simple statistics are first 
computed, such as mean, variance, and the relative frequency distribution, i.e., the histogram. Following this, spectral analysis is 
performed: Individual spectral components are automatically determined, together with their center frequency and associated power, 
i.e., area. Finally, the autospectrum is computed and plotted. Inset of left lower panel: F = frequency in cycles/ beat; P - power in 
msec 2 ; %P = fractional power. In this, as well as in the other R-R interval variability autospectra, the power spectral density (PSD) 
units should be multiplied by JO 3 . 



during tilting of appropriate length were selected for 
analysis. As schematically represented in Figure 1, the 
computer program first calculates the interval tacho- 
gram, i.e., the series of N consecutive R-R intervals, 
and saves them in memory. From sections of tacho- 
gram of 512 interval values, simple statistics (mean 
and variance) of the data are computed. This length of 
the tachogram has been selected as a best compromise 
between the need for a large time series, in order to 
achieve greater accuracy in the computation, and the 
need to obtain stationary recordings, which would be 
easier for short time periods. As an example, for an 
average heart rate of about 70/min, 512 successive 
R-R intervals would amount to approximately 7.5 
minutes. The computer program automatically calcu- 
lates the autoregressive coefficients necessary to de- 
fine the power spectral density estimate (see the Ap- 
pendix). An important feature of the program is that it 
also calculates the model that provides the best statisti- 
cal estimate and prints out the power and frequency of 



every spectral component. Each spectral component is 
presented in absolute units, as well as in normalized 
form, by dividing it by the total power less the DC 
component, if present. This component can also be 
recognized in the graphs of the autospectra on the 
decaying part of the curve near the origin of the abscis- 
sa. Only components greater than 5% of the total pow- 
er were considered significant. Stationarity was as- 
sessed either by pole diagram analysis (see the 
Appendix) or by verifying that a difference of less than 
5% was present between autospectral components cal- 
culated in the two successive 256 beat series constitut- 
ing the whole series of 512 beats. 

In this study, the duration of the periodical phenom- 
ena in the variability signal was measured as a function 
of cardiac beats, rather than seconds. As an example, a 
four-beat periodical component is represented with a 
frequency of 1 :4, i.e., 0.25 cycle/beat. However, this 
frequency is easily converted into hertz equivalents 
(Hz Eq) by dividing it by the average R-R interval 
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length. For instance, if the average R-R length were 
1000 msec, this would correspond to 0.25 Hz Eq. In 
the figures, both units are indicated, whereas, in the 
text, only hertz equivalents are used. 

De Boer et al 23 have compared different methods of 
spectral analysis of the heart rate variability signal and 
found that the analysis of the interval series that repre- 
sents heart period as a function of the beat number not 
only provides comparable results to methods that rep- 
resent heart period as a function of time, but is even 
better suited to the study of the relationship between 
heart period and arterial pressure on a beat-by-beat 
basis. 

After synchronous acquisition and appropriate cali- 
bration, a similar procedure is used to compute the 
spectrum of the arterial pressure data for both systolic 
and diastolic values. 20 

Spectral analysis of the impedance respiratory sig- 
nal was also performed to assess the effects of metro- 
nome breathing on respiratory waveform. 

Statistics 

Data are presented as means ± SE. Student's t test 
was used to determine the significance of the differ- 
ences between rest and tilt. 

Analysis of variance with Scheffe test was used to 
assess the effects of age and of ^-adrenergic receptor 
blockade, as well as the effects of controlled breathing. 

Differences were considered significant at/? < 0.05. 

Regression analysis was used to test the effects of 
age on R-R interval duration variance. 24 



Results 

Age Dependency of Heart Rate Variability 

R-R interval duration variance demonstrated a sig- 
nificant age dependency as it decreased with increas- 
ing age both at rest and during passive upright tilt (90°) 
(Figure 2). This effect appeared significantly repre- 
sented by a curvilinear, i.e., exponential, relationship 
(r = 0.70, p<0.001 at rest, and r = 0.60, /?<0.001 



during tilt). Therefore, subjects are subdivided into 
three age groups of, respectively, 20-30 (n = 30), 
30-45 (n = 10), and 45-60 (n = 17) years (Table 1), 
Furthermore, to account for the possible large differ- 
ences in total power of individual autospectra, vari- 
ance data are presented in absolute units, whereas 
spectral components are normalized by dividing each 
component by total power (less the DC component, if 
present). 

Table 1 indicates that, at rest, there was a slight 
tendency for mean R-R interval to increase with age, 
while R-R variance in the young age group was sig- 
nificantly greater (4097 ± 361 msec 2 ) than both in the 
mid (2581 ±356 msec 2 ) or in the old age group 
(1354 ±205 msec 2 ). 

Effects of Tilt on Heart Rate Variability 

Figure 3 depicts a representative example, in a 
young subject, of the effects of tilting on the time 
series of R-R intervals, i.e., the tachogram, as well as 
on the computed autospectrum. The small fluctuations 
of the instantaneous R-R values around the mean that 
were present both at rest and during tilt. when analyzed 
in their nonrandom components, provided two drasti- 
cally different autospectra. At rest, there were two 
major spectral components at low (—0.1 Hz Eq) and 
high (—0.25 Hz Eq) frequency. The normalized area 
of the low-frequency component (calculated automati- 
cally, see "Materials and Methods" and the Appendix) 
was slightly predominant (58 ± 3%) with an average 
low frequency: high frequency (LF:HF) ratio of 
3.6 ±0.7 (Table 1). Notice that only about 85% of 
total variability is represented by the sum of the LF and 
HF components (Table 1) since in the individual sub- 
jects, smaller components could also be present. 

During tilt, the majority of R-R variability was rep- 
resented by a largely predominant LF component 
(90 ±1%). However, an HF component was also pres- 
ent (9 ± 1%) although at times very small (Figure 3), 
with an average LF:HF ratio of 21 ± 4. As shown in 
Table 1 , heart period during tilting underwent a signifi- 
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Figure 2. Relationship between 
R-R interval variability, expressed 
as variance (o 2 ), and age in the 
study population at rest and during 
90° upright passive tilt. Notice that 
in both cases there is a significant 
exponential relationship (p < 0.001 , 
n = 48). 
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Figure 3. /?-/? interval series, 
i.e., tachogram at rest and during 
passive upright 9Cf tilt. On the auto- 
spectra (bottom panels), two clearly 
separated low- and high-frequency 
components are present at rest. Dur- 
ing tilt, the low-frequency compo- 
nent becomes preponderant. 



0.83 



cant reduction, as expected. However, there was no 
significant correlation between the increases in LF 
component and in heart rate induced by tilting 
(r = 0.l9). 

Autospectra of R-R variability in individual sub- 
jects remained remarkably constant when repeated 
over time. As shown in Table 2, there was no differ- 
ence (/7>0.05), at rest and during tilt, in the results 
obtained in two or three different recording sessions 
with an interval up to 1 year. 

Age Dependency of the Response to Tilt 

During tilt in the old age group, mean R-R was 
significantly longer, and variance smaller, than both in 
the mid and young age groups (Table 1). 



The effects of tilt on the autospectra of R-R interval 
variability when examined in the three different age 
groups demonstrated only minor differences, as LF 
component was smaller and HF component slightly 
greater in the old age group (Table 1). This difference, 
however, was no longer apparent in the LF:HF ratio, 
which was not significantly different in the three age 
groups both at rest and during tilt (Table 1). 

Effects of Controlled Respiration 

The effects of controlled respiration were examined 
in 16 young subjects that were first studied while 
breathing spontaneously (Figure 4, top panels) and, 
subsequently, while breathing following a metronome 
(Figure 4, bottom panels). Respiratory frequency, in 



Table 2, Reproducibility of the Effects of Tilt on R-R Variability 



R-R 
interval 
(msec) 



R-R 
variance 
(msec 2 ) 



Low-frequency component 

Normalized Frequency 
power (Hz Eq) 



High-frequency component 

Normalized Frequency 
power (Hz Eq) 



Rest 

First study (n = 10) 885±43 3264 ±464 62.0±5.0 0.11+0.01 28.5±3.3 0.26±0.02 

Second study (n = 10) 858±32 3508±766 61.4 + 4.6 0.10 + 0.01 29.0±3.7 0.25±0.02 

Third study (n = 4) 914 ±34 4201 +781 56.1±2.6 0.12±0.01 31.3±3.8 0.31^:0.03 

Tilt 

First study (/I = 10) 672 + 24* 2612 + 528 89.5±1.4* 0.09±0.01 6.0 + 0.5* 0.30 ±0.04 

Second study {n — 10) 676 + 33* 1955±303 89.2 + 2.2* 0.10±0.01 6.9±0.07* 0.26±0.03 

Third study (n = 4) 700 ±25* 3381 ±250 87.0 + 2.0* O.IO±0.01 7.7+1.3* 0.27+0.02 

Average delay in time between the first and second study: 138 days (range 7-365). Average delay in time between the second and third 
study: 127 days (range 90-150). 

♦Value during till significantly different from value at rest (p < 0.05). 
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the latter condition, was set at 20/min, i.e., 0.33 Hz, in 
order to separate better the low-frequency from the 
high-frequency respiratory component. 
. This maneuver modified the respiratory waveform, 
as demonstrated by spectral analysis: during spontane- 
ous respiration, nonstationarities in the signal are re- 
flected by a diffuse profile of the spectrum, whereas 
during metronome breathing, the near sinusoidal shape 
of the respiratory waveform is reflected by a single 
narrow component in the spectrum. In either condi- 
tion, the HF component of the spectrum of R-R vari- 
ability coincided with the main respiratory frequency. 
The LF component never coincided with a major com- 
ponent of the spectrum of the respiratory waveform. 



c/b 



J Hz eq 




Figure 4. Representative example of 
autospectra of R-R interval variability 
during spontaneous breathing (top pan- 
els) and during controlled respiration, 
at 20/min (bottom panels). Note that, 
with controlled respiration at rest, the 
relative power of the high-frequency 
component becomes predominant. 



Hz eq 



0.83 



During controlled respiration, at rest, there was a 
significant reduction of the LF component with a con- 
comitant increase in the HF component (Figure 4, bot- 
tom; Table 3). The LF:HF ratio was significantly re- 
duced from 2.5 ±0.3 to 0.7 ±0.1. As the subjects 
were allowed to adjust their tidal volume, during spon- 
taneous respiration the tidal volume (n = 7) was 
552 ± 98 ml, and slightly less (533 ± 76 ml) during 
controlled breathing at 20/min. During metronome 
breathing, there was only a small but not significant 
change in transcutaneous Pco 2 and Po 2 (respectively, 
42 ± 2, 83 ± 4 mm Hg during spontaneous and 
38 ± 2, 79 ± 6 mm Hg during controlled breathing, 
/?>0.05). During tilt, the HF component was de- 



Table 3. Effects of Metronome Breathing on R-R Variability (n ~ J 6) 





R-R 


R-R 


Low-frequency component 


High-frequency component 




interval 
(msec) 


variance 
(msec 2 ) 


Normalized 
power 


Frequency 
(Hz Eq) 


Normalized 
power 


Frequency 
(Hz Eq) 


Rest 














Free resp frequency Hz 0.26 ± 0.02 


884 ±27 


4095 ±436 


56.2±3.9 
* 


0.11 ±0.01 


27.3±2.7 


0.24 ±0.02 


Control resp frequency Hz 0.33 
Tilt 


934 ±29 


4220 ±541 


28.9±4.3 


0.11 ±0.01 


51.6±4.4 


* 

0,33±0.01 


Free resp frequency Hz 0.26 ± 0.03 


678±21t 


2647 ± 308 t 


86.4±2.9t 
* 


0.10±0.01 


7.8±1.4t 


0.25 ±0.02 


Control resp frequency Hz 0.33 


705±22t 


2U9±428t 


71.6±5.3t 


0.10±0.0l 


* 

15.3±3.0t 


0.31 ±0.02 


•Significant difference (p < 0.05). 

tValue during tilt significantly different from value at rest (p < 0.05). 



Pagani et a\ Heart Rate and Arterial Pressure Variabilities 



185 



creased while the LF component was increased, as 
expected; however, this LF increase was less pro- 
nounced than during spontaneous breathing (Figure 4). 
Consequently, the increase in LF:HF ratio induced by 
tilt was smaller during controlled respiration 
(8.3 ±1.6) than during spontaneous breathing 
(16.0 ±3.0). 

Furthermore, during tilt, tidal volume was slightly 
smaller with controlled respiration than with spontane- 
ous breathing (574 ± 100 and 657 ± 48 ml, respective- 
ly), and no significant change (p > 0.05) was observed 
in transcutaneous Pco 2 or Po 2 (respectively, 39 ± 4, 
87 ± 5 mm Hg with spontaneous and 43 ± 3, 82 ± 4 
mm Hg with metronome breathing). 

/3-Adrenergic Receptor Blockade 

The effects of ^-adrenergic receptor blockade were 
assessed in two sets of experiments. Acute /3-adrener- 
gic receptor blockade was obtained in 10 young sub- 
jects by IV bolus of 0.2 mg/kg of propranolol that 
reduced significantly heart rate and arterial blood pres- 
sure (Table 1). Under these conditions, R-R interval 
variance was significantly increased from control (Fig- 
ure 5; Table 1), but normalized autospectral compo- 
nents were not modified significantly (Table 1). Dur- 



ing tilt, there was a significantly lesser reduction in 
R-R interval and a smaller increase in the LF compo- 
nent and in the LF:HF ratio (Figure 5; Table 1). 

Chronic /^-adrenergic receptor blockade was ob- 
tained in 12 young subjects with 0.6 mg/kg p.o. pro- 
pranolol t.i.d. for 6 days, which reduced significantly 
resting heart rate and arterial blood pressure. The ef- 
fects on resting R-R interval variability, like acute 
blockade, were characterized by augmented variance 
(Figure 5; Table 1). However, marked changes were 
observed in the resting autospectra: LF was smaller, 
HF greater, and hence, LF:HF significantly reduced, 
compared with controls (Table 1; Figure 5). During 
tilt, the LF component was reduced and HF component 
increased, compared with control. Thus, during chron- 
ic ^-adrenergic receptor blockade, LF:HF ratio in- 
creased to only 3.78 ± 0.93, which was smaller than 
that observed in the young age group without /3-block- 
ade during tilt (20.79 ± 3.68) but was similar to their 
resting value (3.62 ± 0.70). 

Simultaneous R-R and Blood Pressure Variabilities 

In a group of 7 normotensive subjects (20-60 years 
old) in whom systemic arterial pressure was recorded 
with direct high-fidelity techniques (see "Materials 
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Figure 5. Autospectra of R—R interval variability of a young subject at rest (top panels) and during tilt (bottom panels) under 
control conditions (left panels) and during acute (0.2 mg/kg IV propranolol, middle panels), and chronic (0.6 mg/kg p.o. propranolol 
t.i.d.) fi-adrenergic receptor blockade. Note the progressive reduction in the relative power of the low-frequency component during 
tilt. 
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and Methods") simultaneously with the ECG, the anal- 
ysis of the variability of systolic and diastolic beat-by- 
beat values was also performed. 

Systolic arterial pressure variability was character- 
ized by small oscillations in the time series of beat-by- 
beat values (Figure 6), resulting in a tracing compara- 
ble to the interval tachogram. 

Moreover, the autospectra of systolic arterial pres- 
sure variability demonstrated two major components, 
respectively, a LF of 41 ± 8% and HF of 19 ± 5%, 
with a relative ratio of 2.8 ± 0.7 at rest. During tilt, the 
LF component became largely predominant (73 ± 
8%), with an increase in the LF:HF ratio of 17 ± 5. 
Spectral analysis of diastolic blood pressure gave simi- 
lar results (Table 4). It should be appreciated that the 
simultaneous autospectra of R-R variability (Figure 6; 
Table 4) provided results resembling those obtained 
from the analysis of arterial pressure variability, with 
two components recognized at the same frequencies, 
with similar relative ratios. 

As shown in Figure 7, cross -spectral analysis of 
R-R and blood pressure variability confirmed this 
evaluation by showing coherence only between the 
same 1.F and HF components both at rest and during 
tilt. The phase difference between R-R and arterial 
pressure variability was approximately 0° at HF. At 
LF, the phase presented different patterns in rest and in 
tilt conditions: In the former case, the phase had a 
linear relation with frequency, whereas, in the latter 
case, it oscillated around a fixed value. In both cases, 
the phase corresponding to the central frequency of the 
LF peak was about 60°, i.e. , /s of the entire 360° cycle. 
As LF corresponds to a period of about 10 beats, the 
calculated phase difference amounted to l A of that peri- 
od, i.e., to a delay of 1 .7 beats, with pressure leading. 

R-R and Arterial Pressure Variability in 
Conscious Dogs 

To assess more directly the effects of sympathetic 
innervation on heart rate variability, we performed 
experiments on a group of conscious dogs before and 
after bilateral stellectomy. 

As shown in Figure 8 and Table 5, under control 
conditions R-R variability was almost solely repre- 



sented by a HF respiratory- linked component, whereas 
a very small (8 ± 2%) LF component was present in 
only 50% of the animals. Similarly (Table 5), the auto- 
spectra of both systolic and diastolic blood pressure 
variabilities demonstrated a major HF component. 
During moderate hypotension ( — 9 ± 2% from 85 ± 3 
mm Hg mean arterial pressure) obtained by a continu- 
ous infusion of IV nitroglycerin (32 A-ig/kg per minute), 
heart rate increased 47 ± 9% from 82 ± 5 beats/min as 
a consequence of sympathetic activation. 

Under those conditions, not only was there a signifi- 
cant reduction of total power of R-R variability, but 
also the HF component of the autospectrum was re- 
duced to 42 ± 3%, whereas, in all animals, a LF com- 
ponent of similar power (36 ± 6%) became evident. 
Similar changes were observed in arterial pressure 
autospectra (Table 5). Bilateral stellectomy did not 
significantly modify arterial pressure, heart rate, or 
R— R or pressure variability, as expressed both by vari- 
ance or by their autospectra (Figure 8; Table 5). The IV 
infusion of nitroglycerin, repeated at the same dose, 
reduced arterial pressure and increased heart rate to a 
similar extent (Table 5). The response of R-R variabil- 
ity, however, was modified. Although total power was 
reduced to a value similar to that observed in the intact 
animals, only a HF component could be observed in 
the autospectrum. Conversely, the increase in LF com- 
ponents of pressure autospectra was essentially pre- 
served (Table 5). 

Discussion 

This study in man examines the beat-to-beat oscilla- 
tions which characterize heart rate and arterial pressure 
under various steady state conditions, in the hypothesis 
that the quantitative information provided by the spec- 
tral analysis of these oscillations reflects the interac- 
tion between sympathetic and parasympathetic regula- 
tory activities. 

Various approaches have been proposed to evaluate 
the contributions of parasympathetic and sympathetic 
discharges 25 to the heart rate variability. Parasympa- 
thetic activity has been clinically inferred from peak- 
to-peak variations of the heart period either in the 
clinical laboratory environment 26 or with Holter moni- 



Table 4. Effects of Tilt on Simultaneous R-R Interval and Arterial Pressure Variabilities (n = 7) 









Low-frequency component 


High-frequency component 




Mean value 


Variance 


Normalized 
power 


Frequency 
(Hz Eq) 


Normalized 
power 


Frequency 
(Hz Eq) 


Rest 














R-R 


748 ± 20 msec 


1403 ±571 msec 2 


41.2±8.7 


0.12±0.01 


19.4±5.3 


0.33±0.03 


SAP 


125 ±5 mm Hg 


20 ±3 mm Hg 2 


46.7 ±7.2 


0.11 ±0.02 


27.5±9.0 


0.32±0.03 


DAP 
Tilt 


70 ±5 mm Hg 


15 ±3 mm Hg 2 


48.5 ± 11.7 


0.10±0.01 


31.1 ±7.7 


0.35 ±0.04 


R-R 


633 ±29* msec 


1483 ± 643 msec 2 


73.4 ±8.6* 


0.09±0.01 


7.4±1.9* 


0.33±0.02 


SAP 


129 ±6 mm Hg 


47 ± 10* mm Hg 2 


78.7±4.1* 


0.09±0.01 


9.2±2.1 


0.32±0.03 


DAP 


74 ±4 mm Hg 


29 ±7* mm Hg 2 


73.3±2.5* 


0.08±0.0I 


9.9±1.7+ 


0.35±0.03 



SAP = systolic arterial pressure; DAP = diastolic arterial pressure. 
•Value during tilt significantly different from value at rest (p < 0.05). 
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Figure 6. Simultaneous analysis ofR-R in- 
terval and systolic arterial pressure variabili- 
ties in a young subject. Note the reduction in 
average R-R interval and the slight increase in 
average systolic arterial pressure on passing 
from rest to tilt. The two components of the 
autospectrum of systolic arterial pressure vari- 
ability are similar to those of R-R interval vari- 
ability at rest and during tilt. 
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toring. 27 A different approach has been used by oth- 
ers 1012 ' 5 who utilized spectral techniques to assess the 
frequency components of the heart rate variability sig- 
nal. A possible major advantage of spectral analysis is 
the observation that changes in the sympathetic activ- 
ity to the heart can also be recognized, and hence, 
some index of the instantaneous balance between sym- 
pathetic and vagal activity can be obtained. 15 It should 
be noted that, since heart rate variability signal is a 
pseudorandom phenomenon, previous studies using 
the fast Fourier transform to compute the power spec- 
trum have some technical limitations (see the Appen- 
dix). They include the deterministic nature of the algo- 
rithms used, which, in principle, are applicable only to 
periodical phenomena, 16 the need of windowing the 
data, and the difficulty in defining with certainty the 
relative power of the various spectral components. 

In the present study, some of these limitations were 
avoided because the autoregressive algorithms, which 
are applicable to nonperiodical phenomena, 28 comput- 
ed automatically the number and relative power of the 



various spectral components. Moreover, there was no 
need for windowing or filtering the data. 13 - 22 

Spectral Analysis of Heart Rate Variability and 
Parasympathetic Activity 

In our analysis, we observed at rest two consistent 
major spectral components, a LF at —0.1 Hz and an 
HF at —0.25 Hz. The LF component seems to corre- 
spond to the Mayer waves, 517 while the HF component 
is synchronous with the respiration and has been con- 
sidered as a quantitative evaluation of respiratory ar- 
rhythmia. 13 Since the HF component disappears after 
atropine, it could represent a clinically useful index of 
vagal activity. 15 

In our study, during spontaneous quiet breathing at 
rest, a relatively small area was found to be associated 
with the HF component. However, Pomeranz et al 15 
found a predominant HF component in 8 subjects who 
breathed following a metronome at 15/min. A similar- 
ly enhanced HF component was found in this study in 
16 subjects who breathed following a metronome at 
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Figure 7. Cross-spectral analysis 
(coherence, thick lines; phase, thin 
lines) between R-R interval and sys- 
tolic arterial pressure variabilities 
computed from the data of Table 4. 
Note that, both at rest and during 
tilt, only two spectral components 
show a coherence greater than 0.5 in 
both conditions. Phase is about 0° at 
the respiratory component and nega- 
tive, i.e., pressure leads R—R inter- 
val, at the 0.1 -Hz component. 
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20/min. The same subjects showed a predominant LF 
component when breathing spontaneously. Changes in 
the amount of respiratory arrhythmia and, hence, of 
the HF component could reflect changes in volume, as 
well as in frequency of respiration. 29 * 30 Important mod- 
ulatory mechanisms could also be initiated by changes 
in arterial Po 2 and Pco 2 , and thereby changes in neural 
efferent activities, 31 as well as by variations in arterial 
pressure and consequent alterations in baroreceptor 
input. 32 

In this study, we could not define the exact mecha- 



nisms responsible for the observed increase in HF 
component with controlled breathing, as tidal volume, 
transcutaneous Po 2 and Pco 2 , and arterial pressure did 
not change significantly. 

However, an important role is likely to be played by 
the marked change in respiratory waveform, 33 acting 
either through a more efficient stimulation of lung re- 
ceptors 34 or through fluctuations induced in arterial 
pressure. 32,35 In either case, the increased HF compo- 
nent of R-R variability would indicate that during met- 
ronome breathing there is enhanced vagal tone. 



Table 5. 



Effects of Moderate H ypotension on Simultaneous RR Interval and Arterial Pressure Variabilities in Conscious Dogs 

Low-frequency component High-frequency component 



Mean value 



Variance 



Normalized 
power 



Frequency 
(Hz Eg) 



Normalized 
power 



Frequency 
(Hz Eq) 



Part A 

Intact (n = 8) control 

R-R 755 ± 62 msec 

SAP 119± 18 mm Hg 

DAP 71 ±3 mm Hg 

Nitroglycerin ■ 

(32 jig/kg per min) 

R-R 521 ±37* msec 

SAP 104 ±5* mm Hg 

DAP 66 ±2 mm Hg 
Part B 

Stellectomy (n = 5) control 

R-R 790 ±41 msec 

SAP 137 ± 10 mm Hg 

DAP 76 ±3 mm Hg 

Nitroglycerin 

(32 {Jig/kg per min) 

R-R 617 ± 19* msec 

SAP 110 ±7* mm Hg 

DAP 73 ±4 mm Hg 



58237± 17754 msec 2 4±2 0.14±0.01 

217±112 mm Hg 2 19±4 0.12±0.02 

161 ±71 mm Hg 2 15±5 0.10±0.03 

7513± 3941* msec 2 36 ±6* 0.12±0.01 

45 ±19 mm Hg 2 55 ±4* 0.09 ±0.01 

50±21 mm Hg 2 59 ±7* 0.10±0.01 

42786 ± 14370 msec 2 .3 ±2 0.10 ±0.01 

67 ±28 mm Hg 2 8±3 0.11±0.01 

78 ±22 mm Hg 2 8±4 0.09 ±0.01 

3491 ± 1696* msec 2 0±0 + immeasurable 

16±4 mm Hg 2 38 ±8* 0.08±0.01 

24 ± 13 mm Hg 44 ±6* 0.09±0.01 



63 ±6 
60±4 

64 ±4 

42 ±4* 
27 ±5* 
22 ±5* 



61 ±5 
72 ±4 
63 ±7 



0.24 ±0.02 
0.23 ±0.03 
0.24 ±0.03 



0.26 ±0.03 
0.27 ±0.06 
0.25 ±0.03 



0.23±0.01 
0.25 ±0.01 
0.23 ±0.02 



76±8t 0.24 ±0.02 
46 ±10* 0.26 ±0.02 
37 ±7* 0.25 ±0.02 



SAP - systolic arterial pressure; DAP - diastolic arterial pressure. 

* Value during nitroglycerin infusion significantly different from control (p 

tSignificant difference (p < 0.05) between value in intact and denervated 



< 0.05). 
animals. 
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Spectral Analysis of Heart Rate Variability and 
Sympathetic Activity 

The possibility that spectral analysis could provide 
an index of sympathetic activity is less well estab- 
lished. In animal studies, sympathetic activity has 
been considered either to play no role 10 or else to be 
instrumental in the genesis of the LF rhythm. 12 36 Fur- 
thermore, electroneurographic recordings in dogs 25 
suggest that sympathetic activity might also modulate 
respiratory arrhythmias. This apparent discrepancy of 
results is probably a consequence of the differences in 
the preparations used, i.e., anesthetized and acutely 
decerebrated cats 10 as opposed to conscious dogs. 12 

In our study, we tested the hypothesis that spectral 
analysis of heart rate variability could provide an as- 
sessment of sympathetic tone, by planning experi- 
ments in man in which sympathetic activity was either 
increased functionally or blocked pharmacologically. 
Enhanced sympathetic drive to the heart, as obtained 
by an orthostatic stimulus, was constantly associated 
with a marked increase in the LF and with a decrease in 
the HF component of the autospectrum. We selected a 
passive change of posture in order to minimize the 
influence on heart rate exerted by the muscular effort 
of active standing. 37 

Acute /3-adrenergic receptor blockade with intrave- 
nous propranolol, besides reducing heart rate, had a 
marked effect on R-R variability. At rest variance 
increased, while the relative contribution of LF and HF 
components remained relatively unchanged. During 
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Figure 8. Autospectra of R-R inter- 
val variability of a conscious healthy 
dog in control conditions (left panels) 
and during an IV infusion of nitroglycer- 
in to excite the sympathetic outflow. The 
top panels were obtained with cardiac 
nerves intact, the bottom panels after 
full recovery from bilateral stellectomy. 
Note the presence of a large low-fre- 
quency component during nitroglycerin 
infusion only in the neural ly intact situa- 
tion. 
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tilt, there was a significant reduction of the increase of 
the LF component and of the LF:HF ratio. Similar 
blunting effects exerted by acute /3-adrenergic receptor 
blockade were observed during active standing by 
Pomeranz et al. 15 

After chronic ^-adrenergic receptor blockade, there 
was, at rest, not only an increase in R-R interval vari- 
ance similar to that observed with acute blockade, but 
also a significant effect on spectral components. Thus, 
LF and LF:HF were significantly smaller, and HF was 
significantly greater than under control conditions. 
During tilt, the increase in LF component and in the 
LF:HF ratio was markedly reduced; indeed, the nu- 
merical value of the ratio was similar to that observed 
under resting control conditions without ^-blockade. 

It might be pertinent to recall that Wallin et al 38 
demonstrated, with electroneurographic techniques in 
hypertensive subjects, that chronic /3-adrenergic 
blockade was effective in reducing sympathetic effer- 
ent activity, whereas acute blockade was ineffective. 

As to the effects of metronome breathing on LF, this 
component was reduced at rest, and its increase during 
tilt was blunted, possibly as a consequence of the in- 
hibitory effect on sympathetic rhythms produced by 
pulmonary afferent activity. 34 

As to the effects of aging, it was progressively asso- 
ciated with a reduced R-R variability as expressed by 
variance, and during tilt with a smaller activation of LF 
and a smaller reduction of HF. However, the changes 
in LF:HF ratio during tilt were not significantly differ- 
ent among the various age groups. This approach 
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would suggest that aging is characterized by a new 
equilibrium between the two sections of the visceral 
nervous system rather than by alterations limited to the 
sympathetic outflow. 3940 
Simultaneous Variability of R-R Interval and 
Arterial Pressure 

An additional major observation of this study is that 
the autospectra of R-R interval and of both systolic 
and diastolic arterial pressure beat-to-beat variabilities 
were similar at rest and showed parallel changes with 
tilt. The presence in arterial pressure recordings of an 
HF respiratory component has been traditionally inter- 
preted as a mechanical consquence of respiration, 
which could act directly on intrathoracic vessels or 
indirectly through changes in stroke volume 7,41 and 
heart period. 36 Sympathetic modulation of arterial 
smooth muscle is probably too slow to follow the 0.25- 
Hz respiratory frequency. 42 Obviously, these beat-to- 
beat pressure changes could affect R-R interval 
through complex reflex adjustments, 43 among which 
baroreflexes could have a paramount importance. 44 

The LF components correspond to the well-known 
Mayer 3 waves, a phenomenon that, although described 
in quite artificial experimental conditions, seems to 
pertain to normal human subjects 8 as well. Various 
theories, including myogenic oscillations, central 
rhythms, feedback mechanisms, and "resonance" dis- 
turbances, have been advanced for their interpreta- 
tion. 17 

As to the effects of tilt, arterial pressure beat-by-beat 
variability, as expressed by variance, increased sig- 
nificantly for both systolic and diastolic values. Fur- 
thermore, the LF oscillatory components increased 
signficantly and to an extent similar to that observed in 
the autospectrum of R-R interval, whereas HF compo- 
nents were, in both autospectra, markedly reduced. 

In this respect, it should be mentioned that an in- 
crease in arterial pressure variability has been de- 
scribed also. in essential hypertension, 45 and a greater 
LF component has been documented in daytime as 
opposed to nighttime recordings in ambulatory pa- 
tients. 46 Both conditions are possibly characterized by 
a higher sympathetic activity. 

Cross-spectral analysis of systolic arterial pressure 
and R-R interval variabilities indicated that a high 
degree of coherence existed between the fluctuations 
of these two variables both in recumbency and during 
tilt. In correspondence to the HF component, arterial 
pressure and R-R interval changes occurred in phase, 
whereas each LF pressure change preceded R-R inter- 
val oscillation by about two beats. 

Although this cross-spectral analysis provides no 
direct insight into the mechanisms linking heart period 
and arterial pressure oscillations, with their possible 
neural and non-neural components, 12,47 it supports the 
conclusion that similar information on oscillatory 
rhythms can be obtained from both invasive and nonin- 
vasive studies, not only at rest, but also during aug- 
mented sympathetic activity. 

In a previous study by De Boer et al, 47 only a small 
HF component was observed in the autospectrum of 



diastolic arterial pressure measured in sitting subjects 
either with direct or indirect plethysmography clamp 
method. 48 Technical differences in arterial pressure re- 
cording and the intermediate level of gravitational 
stimulation induced by sitting should account for the 
difference in our data. 

Conscious Dogs 

More direct information on the role of cardiac sym- 
pathetic nerves was derived from a group, of experi- 
ments in conscious dogs before and after bilateral stel- 
lectomy. In the conscious dog, the LF components are, 
under control conditions, either absent 33 or variable 12 - 36 
and small (Table 5). An excitation of sympathetic ac- 
tivity as induced by moderate hypotension obtained 
with phentolamine 36 or, as in our experiments, with 
nitroglycerin, was associated with an increase in LF 
components. This increase was no longer present after 
bilateral stellectomy, which, on the other hand, did not 
modify the LF and HF components characterizing the 
arterial pressure variability (Table 5). Hence, in con- 
scious dogs that had fully recovered from surgery, it 
was possible for the first time to dissociate the LF 
components present in the heart rate and arterial pres- 
sure autospectra by selectively interrupting the cardiac 
sympathetic loop. 

Changes in LF:HF Ratio as Markers of Changes in 
Sympatho-Vagal Balance 

It has been a rather simple but traditional working 
hypothesis that conditions which increase sympathetic 
activity decrease vagal tone and vice versa. 48 How- 
ever, in electrophysiological experiments, recordings 
of the sympathetic and vagal discharges, directed to 
the heart, indicated both a reciprocal 30 and a nonreci- 
procal 31 organization. 

The actual balance between these two outflows is 
likely to have paramount importance not only under 
such physiological conditions as gravitational stimuli, 
exercise, and changes in central command, 51 but also 
during disturbances such as arrhythmias, 52 myocardial 
ischemia, arterial hypertension, 53 or alterations in the 
renin-angiotensin system. 12 

It has recently been suggested that spectral analysis 
of R-R interval variability might reflect this balance. 15 
In our study, the LF: HF ratio appeared to be a conven- 
ient index of such interaction. Indeed, the simulta- 
neous increase in vagal and decrease in sympathetic 
mechanisms produced by metronome breathing were 
reflected by a reduction of LF : HF ratio at rest and by a 
blunting of the increase induced by tilt. Furthermore, 
LF:HF ratio appeared to follow the reductions of sym- 
pathetic activity produced by acute and chronic 
adrenergic receptor blockade both at rest and during 
tilt. Clinical research into those states characterized by 
a disturbance in the neural regulatory mechanisms will 
probably revea! in the near future whether this ap- 
proach will have true pathophysiological significance. 

Appendix 

The basic assumption underlying the proposed signal pro- 
cessing methods is that heart rate and systolic and diastolic 
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blood pressure values fluctuate cycle by cycle, even in stable 
conditions, around a given mean value. Beat-to-beat heart rate 
and arterial blood pressure variability signals are obtained first 
from the original biological signals as a discrete time series y(k) 
where k is the progressive number of detected events. 

Variability signals are intrinsically pseudo-random and can 
be considered as the realization of a stochastic process y(k) 
which, is the output of a linear time-invariant system driven by 
white noise w(k), which constitutes the random component of 
the model and is fully characterized by a mean value which is 
zero for simplicity of notation and a variance X 2 [i.e., 
w(k) = WN(0,X 2 )] according to the model 



w(k)- 



H(z) 



y(k). 



where H(z) is the transfer function in the complex variable z, 
which is determined in a quantitative way and is the model of the 
signal-generating mechanism. A simple structure of input/out- 
put relation is given by the autoregressive (AR) modelization 

y(k)= £ a(i)y(k - i) + w(k) 

i = 1 

where a(i) are the p unknown parameters of indent ifi cation. 
More complex models are described by Box and Jenkins. 28 
Obviously, 



H(z) = 1 - 2 a(i)-z~ 



which relates the transfer function with the identification coeffi- 
cients. Hence, given N samples of variability signals, the prob- 
lem is to estimate H(z) or, more precisely, the vector 8 of 

parameters 

B - | a(l), a<2), a(3) ..... a(p), X 2 \ 

The algorithm applied for such an identification is the Batch 
least squares method via Lcvinson-Durbin recursion. 18 

Two tests are used to check the validity of the assumed 
model. The first one is Anderson's test, which measures the 
whiteness of the prediction error: The identification is not ac- 
cepted if the test is not satisfied within 5% confidence interval. 
After fulfilling Anderson's test, the order of the model is chosen 
as the one which minimizes Akaike's final prediction error 
(FPE) figure of merit. 54 In this way, the model is completely 
determined by order p and vector 8. 

Power spectral density (PSD) estimation P(f) is obtained from 
the following relation 



P(0 = M 



1 - 



p 
1 



a(i)exp(— j27rfi) 



where f is the frequency and the sampling period is unitary. The 
resulting PSD satisfies the criterion of maximum entropy 16 and 
presents many advantages in respect to the methods based on 
classical Fourier analysis (FFT algorithms), 18 namely; a more 
consistent and smoother spectral estimation, a spectral resolu- 
tion which is independent of the number N of samples, and the 
possibility of avoiding windowing procedures. 

Another important advantage of AR modeling with respect to 
the more traditional techniques is the possibility of decompos- 
ing the autocorrelation function and, hence, its transform in the 
frequency domain (i.e., the power spectrum density) in single 
components by applying the residuals theorem. 55 In this way, 
the AR spectrum also provides the individual spectral compo- 
nents in terms of center frequency and of the corresponding 
power in absolute, fractional (see Figure 1), and normalized 
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values. The spectrum is calculated on 512 consecutive cardiac 
beats; stationarity has been tested in two ways: (1) Calculation 
of mean values and variance of variability signals on consecu- 
tive records, by verifying that the values are inside a confidence 
level of 5%, (2) construction of the pole diagram of the auto- 
regressive model, together with the confidence interval in the 
pole estimation, by verifying that consecutive records are char- 
acterized by poles which are inside the fitted confidence interval 
(wide-sense stationarity). The spectra which are presented in 
this paper satisfy both of the preceding conditions. 

With the AR spectral estimation presented in this paper, it is 
possible to improve the statistical consistency, even of very low 
frequency components (less than 0.02 Hz Eq) which correspond 
to oscillations in the signal having a period comparable or supe- 
rior to the considered number of cardiac beats (512). These 
rhythms are always present in the signals and may have an 
important physiological meaning, the analysis of which is out- 
side the aim of the present paper. The method of variability 
signal processing presented has the advantage of decreasing the 
estimation variance of their power. 

The multivariate analysis (cross -spectra and phase spectra) 
between R-R interval and systemic arterial pressure variabili- 
ties are obtained through the calculation of the complex cross- 
spectrum P Ry (f) where x(k) and y(k) are two synchronous dis- 
crete time series. 

-P X y(0 = X(f) Y*(f)/N - G^ffiexpU^f)] 

where X(f) and Y(f) are the discrete Fourier transforms of the 
series (* denotes the complex conjugate), G xy is the amplitude 
cross-spectrum, and 4> Ky is the phase spectrum. The squared 
coherence k 2 is then obtained as 

k 2 xy (o = G Ry (n 2 /rp x (f)P y (f)]- 

The cross-spectrum is calculated via FFT algorithm, using a 
triangular window (Bartlett method) 56 on successive over- 
lapping records of 64 samples each. Results are averaged over 
the whole set of 512 samples. The applied segmentation is a 
compromise between frequency resolution and estimated con- 
sistency. 57 

Squared coherence has values between 0 and 1 and should be 
considered as analogue to the squared correlation coefficient r 2 
in linear regression analysis. Values over 0.5 indicate a signifi- 
cant phase link between R-R interval and pressure variability 
signals. Thus, only spectral components with high coherence 
demonstrate a stable phase-shift between instantaneous R-R 
interval and arterial blood pressure. Its value can be evaluated as 
a function of frequency. 
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